The Reionization of Carbon

Observations suggest that CII was more abundant than CIV in the intergalactic medium towards the end of the hydrogen reionization epoch. This transition provides a unique opportunity to study the enrichment history of intergalactic gas and the growth of the ionizing background (UVB) at early times. We study how carbon absorption evolves from z=10-5 using a cosmological hydrodynamic simulation that includes a self-consistent multifrequency UVB as well as a well-constrained model for galactic outflows to disperse metals. Our predicted UVB is within 2-4 times that of Haardt & Madau (2012), which is fair agreement given the uncertainties. Nonetheless, we use a calibration in post-processing to account for Lyman-alpha forest measurements while preserving the predicted spectral slope and inhomogeneity. The UVB fluctuates spatially in such a way that it always exceeds the volume average in regions where metals are found. This implies both that a spatially-uniform UVB is a poor approximation and that metal absorption is not sensitive to the epoch when HII regions overlap globally even at column densites of 10^{12} cm^{-2}. We find, consistent with observations, that the CII mass fraction drops to low redshift while CIV rises owing the combined effects of a growing UVB and continued addition of carbon in low-density regions. This is mimicked in absorption statistics, which broadly agree with observations at z=6-3 while predicting that the absorber column density distributions rise steeply to the lowest observable columns. Our model reproduces the large observed scatter in the number of low-ionization absorbers per sightline, implying that the scatter does not indicate a partially-neutral Universe at z=6.Comment: 16 pages, 14 figures, accepted to MNRA

Direct Detection of Galactic Halo Dark Matter

The Milky Way Galaxy contains a large, spherical component which is believed to harbor a substantial amount of unseen matter. Recent observations indirectly suggest that as much as half of this ``dark matter'' may be in the form of old, very cool white dwarfs, the remnants of an ancient population of stars as old as the Galaxy itself. We conducted a survey to find faint, cool white dwarfs with large space velocities, indicative of their membership in the Galaxy's spherical halo component. The survey reveals a substantial, directly observed population of old white dwarfs, too faint to be seen in previous surveys. This newly discovered population accounts for at least 2% of the halo dark matter. It provides a natural explanation for the indirect observations, and represents a direct detection of Galactic halo dark matter.Comment: 13 pages, 4 figures, 1 table. Note added after Science Express online publication: This text reflects the correction of a few typographical errors in the online version of the table. It also includes the new constraint on the calculation of d_max which accounts for the fact that the survey could not have detected stars with proper motions below 0.33 arcseconds per year. Published online at ScienceExpress www.sciencemag.org 22 March 2001; 10.1126/science.1059954; To appear in Science 27 April 200

Non-equilibrium chemistry and cooling in the diffuse interstellar medium - II. Shielded gas

We extend the non-equilibrium model for the chemical and thermal evolution of diffuse interstellar gas presented in Richings et al. to account for shielding from the UV radiation field. We attenuate the photochemical rates by dust and by gas, including absorption by HI, H2, HeI, HeII and CO where appropriate. We then use this model to investigate the dominant cooling and heating processes in interstellar gas as it becomes shielded from the UV radiation. We consider a one-dimensional plane-parallel slab of gas irradiated by the interstellar radiation field, either at constant density and temperature or in thermal and pressure equilibrium. The dominant thermal processes tend to form three distinct regions in the clouds. At low column densities, cooling is dominated by ionized metals such as Si II, FeII, FeIII and C II, which are balanced by photoheating, primarily from HI. Once the hydrogen-ionizing radiation becomes attenuated by neutral hydrogen, photoelectric dust heating dominates, while C II becomes dominant for cooling. Finally, dust shielding triggers the formation of CO and suppresses photoelectric heating. The dominant coolants in this fully shielded region are H2 and CO. The column density of the HI-H2 transition predicted by our model is lower at higher density (or at higher pressure for gas clouds in pressure equilibrium) and at higher metallicity, in agreement with previous photodissociation region models. We also compare the HI-H2 transition in our model to two prescriptions for molecular hydrogen formation that have been implemented in hydrodynamic simulations

The Structure of High Strehl Ratio Point-Spread Functions

We describe the symmetries present in the point-spread function (PSF) of an optical system either located in space or corrected by an adaptive o to Strehl ratios of about 70% and higher. We present a formalism for expanding the PSF to arbitrary order in terms of powers of the Fourier transform of the residual phase error, over an arbitrarily shaped and apodized entrance aperture. For traditional unapodized apertures at high Strehl ratios, bright speckles pinned to the bright Airy rings are part of an antisymmetric perturbation of the perfect PSF, arising from the term that is first order in the residual phase error. There are two symmetric second degree terms. One is negative at the center, and, like the first order term, is modulated by the perfect image's field strength -- it reduces to the Marechal approximation at the center of the PSF. The other is non-negative everywhere, zero at the image center, and can be responsible for an extended halo -- which limits the dynamic range of faint companion detection in the darkest portions of the image. In regimes where one or the other term dominates the speckles in an image, the symmetry of the dominant term can be exploited to reduce the effect of those speckles, potentially by an order of magnitude or more. We demonstrate the effects of both secondary obscuration and pupil apodization on the structure of residual speckles, and discuss how these symmetries can be exploited by appropriate telescope and instrument design, observing strategies, and filter bandwidths to improve the dynamic range of high dynamic range AO and space-based observations. Finally, we show that our analysis is relevant to high dynamic range coronagraphy.Comment: Accepted for publication in ApJ; 20 pages, 4 figure

Atmospheric, Evolutionary, and Spectral Models of the Brown Dwarf Gliese 229 B

Theoretical spectra and evolutionary models that span the giant planet--brown dwarf continuum have been computed based on the recent discovery of the brown dwarf, Gliese 229 B. A flux enhancement in the 4--5 micron window is a universal feature from Jovian planets to brown dwarfs. We confirm the existence of methane and water in Gl 229 B's spectrum and find its mass to be 30 to 55 Jovian masses. Although these calculations focus on Gliese 229 B, they are also meant to guide future searches for extra-solar giant planets and brown dwarfs.Comment: 8 pages, plain TeX, plus four postscript figures, gzipped and uuencoded, accepted for Scienc

The Simulated HI Sky at low redshift

Observations of intergalactic neutral hydrogen can provide a wealth of information about structure and galaxy formation, potentially tracing accretion and feedback processes on Mpc scales. Below a column density of NHI ~ 10^19 cm-2, the "edge" or typical observational limit for HI emission from galaxies, simulations predict a cosmic web of extended emission and filamentary structures. We study the distribution of neutral hydrogen and its 21cm emission properties in a cosmological hydrodynamic simulation, to gain more insights into the distribution of HI below NHI ~ 10^19 cm-2. Such Lyman Limit systems are expected to trace out the cosmic web, and are relatively unexplored. Beginning with a 32 h^-1 Mpc simulation, we extract the neutral hydrogen component by determining the neutral fraction, including a post-processed correction for self-shielding based on the thermal pressure. We take into account molecular hydrogen, assuming an average density ratio Omega_H2 / Omega_HI = 0.3 at z = 0. The statistical properties of the HI emission are compared with observations, to assess the reliability of the simulation. The simulated HI distribution robustly describes the full column density range between NHI ~ 10^14 and NHI ~ 10^21 cm-2 and agrees very well with available measurements from observations. Furthermore there is good correspondence in the statistics when looking at the two-point correlation function and the HI mass function. The reconstructed maps are used to simulate observations of existing and future telescopes by adding noise and taking account of the sensitivity of the telescopes. The general agreement in statistical properties of HI suggests that neutral hydrogen as modeled in this hydrodynamic simulation is a fair representation of that in the Universe. (abridged)Comment: 20 pages, 17 figures, Accepted for publication in A&A, figures compressed to low resolution; high-resolution version available at: http://www.astro.rug.nl/~popping/simulated_HI_sky.pd

Metal Enrichment in the Reionization Epoch

The presence of elements heavier than helium ("metals") is of fundamental importance for a large number of astrophysical processes occurring in planet, star and galaxy formation; it also affects cosmic structure formation and evolution in several ways. Even a small amount of heavy elements can dramatically alter the chemistry of the gas, opening the path to complex molecules. Metals might enhance the ability of the gas to radiate away its thermal energy, thus favoring the formation of gravitationally bound objects; they can also condensate in a solid phase (dust grains), partly or totally blocking radiation from luminous sources. Finally, they represent useful tracers of energy deposition by stars and probe the physical properties of the environment by absorption or emission lines. Last, but certainly not least, life -- as we know it on Earth -- is tightly related to the presence of at least some of the heavy elements. In this pedagogical review I will concentrate on the connection between early metal enrichment and cosmic reionization. As we will see these two processes are intimately connected and their joint study might turn out to be fundamental in understanding the overall evolution of the Universe during the first billion years after the Big Bang, an epoch corresponding to redshifts z>6.Comment: Book chapter in Understanding the Epoch of Cosmic Reionization: Challenges and Progress, Springer International Publishing, Ed. Andrei Mesinger, ISBN 978-3-319-21956-1. arXiv admin note: text overlap with arXiv:astro-ph/0007248 by other author

The lensing properties of subhaloes in massive elliptical galaxies in sterile neutrino cosmologies

We use high-resolution hydrodynamical simulations run with the EAGLE model of galaxy formation to study the differences between the properties of - and subsequently the lensing signal from - subhaloes of massive elliptical galaxies at redshift 0.2, in Cold and Sterile Neutrino (SN) Dark Matter models. We focus on the two 7 keV SN models that bracket the range of matter power spectra compatible with resonantly produced SN as the source of the observed 3.5 keV line. We derive an accurate parametrization for the subhalo mass function in these two SN models relative to cold dark matter (CDM), as well as the subhalo spatial distribution, density profile, and projected number density and the dark matter fraction in subhaloes. We create mock lensing maps from the simulated haloes to study the differences in the lensing signal in the framework of subhalo detection. We find that subhalo convergence is well described by a lognormal distribution and that signal of subhaloes in the power spectrum is lower in SN models with respect to CDM, at a level of 10-80 per cent, depending on the scale. However, the scatter between different projections is large and might make the use of power spectrum studies on the typical scales of current lensing images very difficult. Moreover, in the framework of individual detections through gravitational imaging a sample of ≃30 lenses with an average sensitivity of Msub = 5 × 107 M☉ would be required to discriminate between CDM and the considered sterile neutrino models